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Measurement of a Sign-Changing Two-Gap Superconducting Phase in Electron-Doped 

Ba(Fe! x Co x ) 2 As2 Single Crystals using Scanning Tunneling Spectroscopy 

M. L. TEAGUE 1 , G. K. DRAYNA 1 , G. P. LOCKHART 1 , P. CHENG 2 , B. SHEN 2 , H.-H. WEN 2 and N.-C. YEH 1 
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2 Institute of Physics, Chinese Academy of Sciences, China 
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Scanning tunneling spectroscopic studies of Ba(Fei_ x Co x )2As2 (x = 0.06, 0.12) single crystals 
reveal direct evidence for predominantly two-gap superconductivity. These gaps decrease with in- 
creasing temperature and vanish above the superconducting transition T c . The two-gap nature 
and the slightly doping- and energy-dependent quasiparticle scattering interferences near the wave- 
vectors (±7r, 0) and (0, ±7r) are consistent with sign-changing s-wave superconductivity. The excess 
zero-bias conductance and the large gap-to- T c ratios suggest dominant unitary impurity scattering. 

PACS numbers: 74.55. +V, 74.70.Xa, 74.25.Jb 



The recent discovery of iron-based superconductors [l|- 
1 101 ] has renewed intense research activities in supercon- 
ductivity. Comparison of the similarities and contrasts 
between the iron-based compounds and the cuprates can 
provide useful insights into the microscopic mechanism 
for high-temperature superconductivity |ll|. In particu- 
lar, antiferromagnetic spin fluctuations appear to influ- 
ence the physical properties of both cuprate and iron- 
based superconductors IllMloTl . On the other hand, the- 
oretical calculations |8j, Il5l4l7| and angle resolved pho- 
toemission spectroscopy (ARPES) [18 20] suggest that 
multi-bands and inter-Fermi surface interactions are cru- 
cial to superconductivity in the iron-based compounds. 

Whil e subst antial experimental results from ARPES 
studies 



18 2 



and phase sensitive measurements of var- 
ious single crystalline iron-based compounds [2l|, ]22j are 
supportive of the scenario of two-gap superconductivity 
with sign-changing s-wave (s ) order parameters for the 
hole and electron Fermi pockets, tunneling and point- 
contact spectroscopic studies of the iron pnictides ap- 
pear to be inconclusive (231428 1. For instance, reports 
of point-contact spectroscopy and scanning tunneling 
spectroscopy (STS) have suggested either BCS-like s- 
wave superconductivity 12311 o r nodes in the supercon- 
ducting order parameter [24J, |28j of the "1111" iron- 
arsenides LnFeAsOi_ x F x (Ln: trivalent rare-earth ele- 
ments). On the other hand, STS studies of the electron 
and hole-doped "122" iron arsenides Ba(Fei_ x Co x )2As2 
and (Bai_ x K x )Fe2As2 have revealed findings ranging 
from strong spatial variations in the tunneling spectra 
with occasional observation of a large superconducting 
gap 25 1 to moderate spatial variations with predomi- 
nantly a small superconducting gap [26|. In the non- 
superconducting limit, nematic surface reconstructions 



have been observed [27] 



In this letter we report direct STS evidence for two- 
gap superconductivity in the electron-doped 122 system 
Ba(Fei_ x Co x )2As2 of two different doping levels. For 
each doping level, two different energy gaps can be clearly 
resolved at T < T c . Both gaps decrease monotonically 



with increasing temperature and then completely vanish 
above T c . The gap values agree favorably with those ob- 
tained from ARPES so that the larger gap Ar may be as- 
sociated with the hole-like Fermi pockets and the smaller 
gap Am with the electron-like Fermi pockets. Moreover, 
Fourier transformation (FT) of the tunneling conduc- 
tance reveals energy and doping dependent quasiparti- 
cle scattering interferences (QPI) near the nesting wave- 
vectors (±7r, 0) and (0, ±tt) between the Fermi pockets 
at r and M, which is consistent with the sign- changing 
order parameters for the hole and electron pockets [l7j . 
Finally, excess zero-bias tunneling conductance and the 
large 2Ar,M/(fcsT c ) ratios for both doping levels may be 
attributed to strong unitary impurity scattering [29J ] . 

The Ba(Fei_ x Co x )2As2 samples investigated in this 
work are single crystals with x = 0.06 (undcrdoped) and 
0.12 (overdoped), and the corresponding superconduct- 
ing transition temperatures are T c = 14 and 20 K, re- 
spectively. The single crystals were grown from the flux 
method [3Pj , and details of the synthesis and characteri- 
zation of the samples have been described elsewhere [30l — 
l32j. Given the reactive nature of Ba(Fei_ x Co x )2As2, 
freshly cleaved surfaces were essential for the STS stud- 
ies. To date all reported STS studies were carried out 
on samples that were mechanically cleaved under ultra- 
high vacuum (UHV) conditions and at cryogenic temper- 
tures [25l427| . and the resulting sample surfaces all ex- 
hibited significant (2 x 1) reconstructions of the Fe(Co)- 
layer. We chose to perform mechanical cleavage of the 
single crystals in pure argon atmosphere at room temper- 
ature, well above the tetragonal-to-orthorhombic struc- 
tural phase transition. The cleaved samples were loaded 
in situ onto the cryogenic probe of our homemade scan- 
ning tunneling microscope (STM) in argon. The sealed 
STM assembly was subsequently evacuated and cooled 
to 6 K in UHV with a base pressure at ~ 10~ 10 Torr. 

Both spatially resolved topography and tunnelling 
conductance (dl/dV) versus energy (E — eV) spec- 
troscopy were acquired pixel-by-pixel over an extended 
area of each sample simultaneously, with tunneling cur- 
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FIG. 1. (color online) Direct spectroscopic evidences for two- 
gap superconductivity in Ba(Fei_ x Co x )2As2: (a) Left panel: 
normalized tunneling conductance (dl /dV) vs. bias voltage 
(V) spectra taken at T — 6, 8 and 15 K for the sample with 
x = 0.06 and T c = 14 K. The solid lines represent theoretical 
fittings to spectra using the Dynes formula in Eq. (1) modified 
for two-gap BCS superconductors. Two distinct tunneling 
gaps Ar and Am can be identified from the spectrum at T 
= 6 K. Right panel: The tunneling gaps Ar and Am as a 
function of the reduced temperature (T/T c ) are shown by the 
symbols and solid lines. The error bars indicate the widths of 
the gap distributions obtained from the fitting using Eq. (1). 
(b) Left panel: (dl/dV) vs. (V) spectra taken at T = 6, 10 
and 21 K for the sample with x = 0.12 and T c = 20 K. Right 
panel: A r ,M-vs.-(T/T c ). 



rents along the crystalline c-axis. The typical junction 
resistance was kept at ~ 1 Gil. To remove slight varia- 
tions in the tunnel junction resistance from pixel to pixel, 
the differential conductance at each pixel is normalized to 
the polynomial fit to its high-energy conductance back- 
ground from \E\ = A max + 1.5 meV to \E\ = A max + 6.5 
meV. Detailed survey of the surface topography and tun- 
neling conductance spectra was carried out over typically 
(5.4 x 5.4) nm 2 and (6.0 x 6.0) nm 2 areas, and each area 
was subdivided into (128 x 128) pixels. Generally the 
tunneling spectra appeared to be relatively consistent 
throughout each scanned area, with representative point 
spectra shown in the left panels of Figs. 1(a) and 1(b) for 
doping levels x = 0.06 and 0.12, respectively. Two pre- 
dominant tunneling gap features are apparent for both 
doping levels. In comparison with the ARPES data [la - 
|20| , we may assign the larger gap to the superconducting 
gap of the hole Fermi surface at the T-point of the Bril- 
louin zone, Ar, and the smaller gap to that of the electron 
Fermi surface at the M-point, Am- However, upon closer 
inspection, we note that the larger gap features often ex- 
hibit broadening or even slight splitting, as exemplified 
in Figs. l(a)-(b) and 3(c). The physical origin of this 
splitting/broadening is unknown. 

Next, we employ a phenomenological fitting general- 



ized from the Dynes formula [33j to analyze the spectra 
and we restrict to two-gap superconductivity in our anal- 
ysis. Specifically, the normalized tunneling conductance 
G for a metal-insulator-superconductor junction in the 
case of a two-gap superconductor may be given by: 



G = A 



E 



i=r,M 



B t / Re 



(E-iTJidf/dE^E-eV 



dE. 



y/(E-iTi)*-£j 

(1) 
Here A and Bi are positive constants, I\ denotes the 
quasiparticle scattering rate associated with the super- 
conducting gap Aj, and f(E) is the Fermi function. 
Hence, by applying Eq. (1) to the temperature dependent 
tunneling conductance in the left panels of Figs. l(a)-(b), 
we obtain temperature dependent values for Ar and Am, 
which are illustrated in the right panels of Figs. 1(a)- 
(b). Both gaps are particle-hole symmetric (see Fig. 2) 
and vanish immediately above T c for both doping levels, 
implying that Ar and Am are indeed superconducting 
gaps. Further, the quasiparticle scattering rates derived 
from Eq. (1) are very large even at T — 6 K, show- 
ing (r r /A r ) = 0.4 and 0.5 for x = 0.06 and 0.12, and 
(r M /A M ) = 0.1 for both x = 0.06 and 0.12. 

It is apparent from Fig. 1 that the two-gap fitting 
is not ideal, which may be attributed to the following. 
First, the generalized Dynes formula does not explicitly 
consider the possibility of a sign-changing s-wave order 
parameter, the latter is theoretically shown to be very 
sensitive to unitary impurities so that the zero-bias con- 
ductance may be strongly enhanced without requiring 
a large (r/A) ratio 29]. Second, there may be different 
gaps associated with the two hole-pockets, so that Eq. (1) 
is not consistent with the detailed electronic structures. 

The results shown in Figs. l(a)-(b) are representa- 
tive of the spectral characteristics of both samples, as 
manifested by the gap spatial maps and the correspond- 
ing histograms for both samples in Figs. 2(a)-(b), where 
the gap values are empirically determined as one-half of 
the peak-to-peak values. Overall (Ar) = 10.0 meV and 
(A M ) = 5.0 meV for x = 0.12 and (A r ) = 8.0 meV and 
(A M ) = 4.0 meV for x = 0.06. 

Our findings of the two-gap spectra differ from pre- 
vious STS studies [25|, [26|, which may be the result 
of differences in the surface preparation. Specifically, 
we show in Fig. 3(a) an example of the atomically re- 
solved surface topography of the overdoped sample over 
a (5.4 x 5.4) nm 2 area. We find that cleaving the samples 
under argon gas at room temperature generally resulted 
in fragmented surfaces, as exemplified by the height his- 
togram in Fig. 3(b) for the same (5.4 x 5.4) nm 2 area. In 
particular, we note that the topography in Fig. 3(a) ex- 
hibited no apparently reconstructed (1x2) surfaces, and 
the overall height variations were limited to within one 
c-axis lattice constant cq = 1.239 nm. This finding is in 
contrast to the mostly flat and reconstructed surfaces re- 
ported by other groups for cold-cleaved samples [251427 1. 
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FIG. 2. (color online) Superconducting gap maps and his- 
tograms at T = 6 K: (a) Left to right: The first two pan- 
els correspond to the Am and Ar maps for the underdoped 
sample (x = 0.06), and the right two panel represent the cor- 
responding histograms for both the quasiparticle (solid bars) 
and quasihole (shaded bars) branches, showing particle-hole 
symmetry and the mean values of (|Am|) = 4 meV and 
(|Ar|) = 8 meV. (b) The left two panels are respectively the 
Am and Ar maps for the sample with x = 0.12. The right two 
panels are histograms of Am and Ar, showing particle- hole 
symmetry and (|Am|) = 5 meV, (|Ar|) = 10 meV. 
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FIG. 3. (color online) Correlation of atomically resolved sur- 
face topography with tunneling conductance spectra: (a) Sur- 
face topography of the overdoped sample (x = 0.12) over a 
(5.4 x 5.4)nm area, (b) The height histogram of the area 
in (a), showing height variations within one lattice constant 
along the c-axis. (c) Spatial evolution of the normalized 
(dl /dV)-vs.-V spectra across a horizontal line slightly below 
the middle of the topography image in (a), which reveals con- 
sistent two-gap features. 



Despite the fragmented surfaces, consistent two-gap spec- 
tral features were found throughout the scanned area, as 
exemplified in Fig. 3(c) for the tunneling spectra taken 
near the middle of Fig. 3(a). 

Next, we examine the spatial variations in the tunnel- 
ing conductance at constant bias voltages. In Figs. 4(a)- 
(b), we illustrate the tunneling conductance maps of two 
samples and under three different bias voltages of V = 0, 
((A M )/e) and «A r )/e). The Fourier transformed (FT) 
tunneling conductance at V = ((A]yr)/e) and ((Ar)/e) 



is shown in Fig. 4(c) for x — 0.06 and in Fig. 4(d) 
for x — 0.12. We find that dominant QPI occur at 
three wave- vectors [22j : q\ between two electron pockets 
across the first Brillouin zone, q\ ~ (±27r, 0)/(0, ±27r); 
qi between the hole- and electron-pockets at T and M- 
points, g2 ~ (± 7r I 0)/(0, ±7r); and 93 between two ad- 
jacent electron pockets, q% ~ (±, 7r, ±7r). For the un- 
derdoped sample, strong QPI at qi and absence of QPI 
at 53 is consistent with the five-orbital theoretical cal- 
culations la, UM, |34| for two possible scenarios: one is 
scalar-impurity QPI between sign-reversing order param- 
eters associated with the hole- and electron-pockets, and 
the other is magnetic impurity QPI between order pa- 
rameters of the same sign. Given that scalar scatterers 
are generally more common than magnetic impurities, we 
suggest that the prevailing QPI wave- vectors at qi are the 
result of sign- reversing order parameters. In this context, 
the appearance of (73 wave- vectors in the overdoped sam- 
ple might be due to a larger density of Co-atoms in the 
surface Fe/Co layers: Unlike those in the bulk, the charge 
transfer from the surface Co-atoms is incomplete so that 
they may behave like magnetic impurities |15lll6ll34ll35j. 

It is worth noting that for both samples the (72 wave- 
vectors appear to exhibit one preferential direction. 
In addition, energy-dependent features at small wave- 
vectors (|g| < 7r/2) are observed along the same direction. 
These findings may be related to the nematic-order wave- 
vectors found in the parent state of these compounds |27| , 
although our limited momentum resolution cannot pro- 
vide detailed comparison. Overall, we attribute q\, 52, 
and qz to the QPI scattering wave-vectors because they 
are not only slightly energy dependent but also doping 
dependent, and so they cannot be simply attributed to 
Bragg diffractions of the lattices. The relatively weak en- 
ergy dependence of the QPI wave- vectors is the result of 
small and nearly isotropic Fermi pockets in the iron ar- 
senides, which differs from the highly energy-dependent 
QPI wave- vectors in the cuprate superconductors [36h39J . 

Finally, we note that the relatively high zero-bias con- 
ductance in all tunneling spectra at T < T c is sugges- 
tive of dominant unitary impurity scattering [291 ]. Fur- 
ther, unitary impurity effects on the suppression of T c 
for the sign-changing s-wave superconductors are found 
to be as significant as those on the d-wave superconduc- 
tors [29(. On the other hand, unitary impurity effects 
on suppressing the sign-changing s-wave order parame- 
ters involve sign-dependent components that are partially 
cancelled, and are therefore weakened J29J. Hence, the 
large ratios of (2(Ar,M))/(fcs7c) for the 122 system, with 
(2(A M ))/(fcsT c ) ~ 6.6 (5.8) and (2<A r » / (k B T c ) ~ 13.2 
(11.6) for x = 0.06 (0.12), may be attributed to signif- 
icant unitary impurity scattering in these sign-changing 
s-wave superconductors. A possible source for the uni- 
tary impurity scattering may be associated with disorder 
in the Co doping into the Fe-planes, which is in contrast 
to the situation in the 1111 iron arsenides where doping 




11.3 



< 
I 0,7 



11.3 



L 



X(nm) 6.0 X(nm) 6.0 

|1.3 I '>;■!. i 



olflHMHB 10.7 
X(nm) 6.0 X(nm) 6.0 




X(nm) 6.0 





FIG. 4. (color online) Tunneling conductance maps of 
Ba(Fei_ x Co x )2As2 samples at three constant bias voltages 

V = 0, ((A M )/e), «A r )/e) and for T = 6 K: (a) x = 0.06 
and (b) x = 0.12. (c) Fourier transformed (FT) tunneling 
conductance maps at V = ((Am)/c) = 4 meV (left panel) and 

V = ((Ar)/e) = 8 meV (right panel) for x = 0.06. (d) FT 
conductance maps at V = ((Aii)/e) = 5 meV (left panel) and 

V = ((A r )/e) = 10 meV (right panel) for x = 0.12. The thin 
white lines define the first and second Brillouin zones, and we 
have used the lattice constant a = 0.53 nm for the length of 
the first Brillouin zone (2ir/a). All three QPI wave- vectors 
51, 52 and 53 are slightly energy and doping dependent. 



takes place in the charge reservior [5J so that the disor- 
der effect is weaker and the (2A)/(fcsT c ) ratio is small 
(~ 3) [23j,|28[, whereas the T c values are generally higher. 
In summary, we have demonstrated direct STS evi- 
dence for two-gap superconductivity in electron-doped 
Ba(Fei_ x Co x )2As2 single crystals of two doping lev- 
els. The Fourier transformed tunneling conductance re- 
veals strong quasiparticle scattering interferences near 
the nesting wave-vector between the hole Fermi pock- 
ets at r and the electron Fermi pockets at M, consistent 
with sign-changing order parameters of the two Fermi 
pockets. The excess zero-bias conductance and the large 
2Ar,M/(fcsT c ) ratios for both doping levels may be at- 
tributed to significant unitary impurity scattering in a 
sign-changing s-wave superconducting system. 
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